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ABSTRACT 

We present a new method for studying the intrinsic brightness temperatures of the parsec-scale jet cores 
of Active Galactic Nuclei (AGN). Our method uses observed superluminal motions and observed brightness 
temperatures for a large sample of AGN to constrain the characteristic intrinsic brightness temperature of the 
sample as a whole. To study changes in intrinsic brightness temperature, we assume that the Doppler factors of 
individual jets are constant in time as justified by their relatively small changes in observed flux density. We find 
that in their median-low brightness temperature state, the sources in our sample have a narrow range of intrinsic 
brightness temperatures centered on a characteristic temperature, 7^,,, ~ 3 x 10'° K, which is close to the value 
expected for equipartition, when the energy in the radiating particles equals the energy stored in the magnetic 
fields. However, in their maximum brightness state, we find that sources in our sample have a characteristic 
intrinsic brightness temperature greater than 2 x 10" K, which is well in excess of the equipartition temperature. 
In this state, we estimate the energy in radiating particles exceeds the energy in the magnetic field by a factor 
of ~ 10^. We suggest that the excess of particle energy when sources are in their maximum brightness state is 
due to injection or acceleration of particles at the base of the jet. Our results suggest that the common method 
of estimating jet Doppler factors by using a single measurement of observed brightness temperature and/or the 
assumption of equipartition may lead to large scatter or systematic errors in the derived values. 

Subject headings: galaxies : active — galaxies: jets — galaxies: kinematics and dynamics — radiation mech- 
anisms: non-thermal — radio continuum: galaxies 



1. INTRODUCTION 

A key physical property of extragalactic radio jets is the re- 
lationship between the energy stored in the radiating particles 
and the energy stored in the magnetic field. Burbidge ( 195^ 
originally suggested that these energies were approximately 
in balance, or equipartition, as a way of minimizing the total 
energy stored in exte nded radio lobes, and more recent work 
iCroston et al J 12005) has indeed shown that extended radio 
lobes are at or very near equipartition. However, it is still un- 
known whether this balance is established at, or close to, the 
base of the jet, or if jets begin far out of equipartition and only 
reach equipartition after significant periods of time. 

The energy balance at the base of extragalactic radio jets 
can be studied by measuring the brightness temperature of 
jet cores using the Very Long Baseline Array (VLBA). The 
intrinsic brightness temperature, T},,,, is related to the energy 
balance via the following expression (R eadhead 1994) : 
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where 77 = Up/uB is the ratio of the energy densities of the 
radiating particles, Up, to the magnetic field, ug- Tgq ~ 5 x 
10'" K is the equipartition brightness temperature defined by 
l^eadhead ( 1994), and the exponent of rj assumes an power- 
law Lorentz factor distribution for the radiating particles with 
index p = 2.5 for N-fd-f = K-f'^'dj. This corresponds to a spec- 
tral index of a = -0.75 {S oc i^^") in the optically thin parts of 
the jet. 

The intrinsic brightness temperatures of AGN cores are lim- 
ited eith er by the inverse-Com pton catastrophe to be Ti„, < 
10'^ K jKellermann & Pau Unv-Toth 1969) or by the possi- 
bilit y that jets are near equ ipartition to begin with , so Tint — 
T,a ("Singal & GoDal- Krishna..l985i iSiDgal 198^, iReadhead 
1994). If AGN cores are instead limited by the inverse- 
Compton catastrophe, th ey must be very far out o f equipar- 
tition with 77 ~ 10" (e.g. lKellermann & Pauliny^oth.1969) . 

Unfortunately, it is very difficult to measure intrinsic bright- 
ness temperatures because AGN jets are highly relativistic 
and therefore Doppler boosted. So with Very Long Baseline 
Interferometry (VLBI), we do not measure the 7^,,, directly, 
but rather the observed brightness temperature'" Tabs = STi,„, 
where S is the Doppler factor of the jet, given by 



l-flcosO 

where f3 is the speed of the relativistic flow in units of the 
speed of light, and 9 is the angle the jet makes with the line 
of sight. The speed of the jet is often characterized by its bulk 
Lorentz factor, T = 1 /^l-/?^. 

For this paper we take all observed brightness temperatures from the 
point of view of an observer co-moving with the AGN host galaxy, so we 
have corrected for redshift. We also assume that radiation is emitted isotrop- 
ically in the co-moving frame, and that the parsec-scale jet cores have a flat 
spectrum {a = 0) due to optical depth effects. Our choice of cosmology is 
Ho = 70 km/s/Mpc, Qm = 0.3, and Ha = 0.7. 
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Fig. 1 . — Plot of 1000 fictional sources from a program which simulates 
a relativistically beamed population of AGN jets selected on the basis of ob- 
served flux. All sources are given the same intrinsic brightness temperature. 
The solid line represents where sources observed at the critical angle would 
lie on this plot. The dashed line represents the possible apparent speeds of a 
r = 30 source with intrinsic brightness temperature given by Ti„, . 

Observed brightness temperatures from ground-based 
VLBI can range up to lower limits of Tghs > 5 x lO''' K (e.g. 
iKovalev et alJ20 05). and observed brightness temperatures in 
excess of 1 x 10'"^ K have been reported with space VLBI 
(e.g. iHoriuchi et al. 2004). Horiuchi et al. (2004) find typ- 
ical observed brightness temperatures in the range of 1 x 
10" K to 1 X 10'^ K from the VSOP 5 GHz survey, and 
IKovalev et al. (2005) find a similar range from ground-based 
VLBA measurements at 15 GHz. The recent discovery of 
very nearby scintillation screens has brought the maximum 
brightness temperatures deduced from intra-day variability 
(IDV) into the same range as those measured through VLBI 
even though interstellar scintillation is potentially sensitive 
to much smaller angular sizes than the VLBA ( Lo ve 11 et al.. 
1^03, and references therein). 

In this paper, we present a new method for deducing the 
typical intrinsic brightness temperature for parsec-scale AGN 
cores from a population of well studied VLBA jets. These 
cores represent the optically thick region at the start of the jet 
as seen in VLBI imag es. We combine obs erved brightness 
temperature data from I Kovalev et all ({ 2005') and proper mo- 
tion data from the 2 cm survev ("e.g.lKellermann et al. 2004)- 
In ^we discuss the theory behind our approach, and in ^ 
we present our results and discuss the implications for energy 
balance between particles and magnetic fields at the base of 
AGN jets. 

2. METHODOLOGY 

The basis of our approach is to compare the observed 
brightness temperatures, Tgbs, for a collection of parsec-scale 
AGN jet cores with the proper motions of their jets, Papp- 
Proper motion, Papp = (3smd/{l - l3cos9), will serve as a 
proxy for the jet Doppler factor, 6, which we cannot easily 
measure. For the purposes of this analysis, we will assume 
that the speeds of the fastest jet components are the same as 
the flow speed through the jet core. We will also assume the 
jets are straight, with no bends between the jet core and the 
location of jet components. To the degree that these assump- 
tions are incorrect, we should expect some scatter in our final 
results. 



To generate a simple model, we make two further assump- 
tions: (1) that all jets in our sample have the same intrinsic 
brightness temperature, 7},,,, and (2) that all sources are close 
to the critical angle, 6^ = arccos f3, for the maximal superlu- 
minal motion at a given /3. Under these strict assumptions, 
S ~ (3„pp and 7;^ ~ (3appTi„t, and therefore I3app is a sim- 
ple function of the observed brightness temperature, as rep- 
resented by the solid line in Figure I . 

Of course, the assumption that all jets are near the critical 
angle is poor. For a flux-limited sample, the majority will 
lie inside the critical angle where the Doppler factor exceeds 
the Lorentz factor ( Vermeulen & Cohen 1994). The probabil- 
ity distribution for jet orientations in a beamed, flux-limited 
sample has been studied in detail b y fVermeulen & Cohe^ 
([1994) and Lister & M arscher ( 1997 ). In Figure 1, (3app is 
plotted versus Tobs from a numerical simulation of a flux- 
limited sample of 1000 AG N jets. This simulation is based on 
iLister & Marsched (IT997I) and assumes that all jets have the 
same intrinsic brightness temperature Tint =5 x 10'" K. The 
jet speeds in the population range from /3 = 0.05 to F = 30, 
and are distributed according to n(T) oc F"'^. These simu- 
lated jets are drawn from a randomly-oriented parent popula- 
tion with an intrinsic (unbeamed) luminosity function consis- 
tent w ith FR-II radio galaxies described in Lister & Marsche3 
([l9^|. The luminosity function has a lower cutoff at 1 x 
10^^ W//z"', and incorporates pure exponential luminosity 
evolution while maintaining a constant co-moving space den- 
sity out to z = 4. 

Approximately 75% of the sources in Figure 1 fall below 
and to the right of the solid line representing the critical angle; 
these sources are inside the critical angle, and while they have 
large Doppler factors, many of them are close enough to the 
line of sight to have small proper motions. The dotted line in 
Figure 1 shows the apparent speed for a jet with the maximal 
Lorentz factor of F = 30 in the simulation as a function of Tubs, 
corresponding to a varying jet viewing angle. The dotted line, 
therefore, represents an upper envelope for the simulation. 

The scatter in Figure 1 gives a realistic idea of what we 
might expect for a "best case" correlation between (3app and 
Tobs for real data under the assumptions that jets are straight, 
pattern speeds = flow speeds, and all jets have the same Ti„,. 
However, even with the scatter induced by selection effects, 
the trend is clear, particularly the lack of sources with both 
low observed brightness temperature and fast proper motion. 

3. RESULTS AND DISCUSSION 

We have collected multi-epoch brightness temperature 
measurements for 106 sources that all have well studied 
proper motions and at least 5 epochs of brightness temper- 
ature measurements. The brightness temperature measure- 
ments are taken from Kovalev et al. (2005), and the proper 
motion measurements are updated from Kellerman n et alJ 
(2004) to include additional epochs from the VLBA" 2cm 
Surve y'^ through 2() 02 (E. Ros et al., in prep.). As described 
in iiKellermann et al] ( 12004 ) we use the fastest proper motion 
for each source, and we only consider those motions that are 
ranked as "Excellent" (E) or "Good" (G) by the criteria laid 
out in that paper Of the 106 sources plotted here, 70 are 
members of the 133 source, flux-limited MOJAVE'^ sample 

' ' The Very Long Baseline Array is operated by the National Radio As- 
tronomy Observatory, a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

http://www.cv.nrao.edu/2cmsurvey/ 

http://www.physics.purdue.edu/astro/MOJAVE/ 
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jLister & HomarJ2005l) . The remaining 35 sources are mem- 
bers of the 2cm Survey, but do not meet the selection criteria 
of the MOJAVE program. 

Figure 2 is a three panel plot of f3app vs Tabs- Figure 2a 
contains 106 fictional sources from the simulation plotted in 
Figure 1. The 106 fictional sources are selected in a manner 
analogous to our actual data, with 70 of them chosen at ran- 
dom from the 133 brightest sources in the simulation and the 
remaining 36 chosen at random from the rest. This panel is 
included to illustrate that selection criteria similar to our real 
sample will still produce a plot very similar to Figure 1 with 
about 75% of the points falling below and to the right of the 
solid line representing sources at the critical angle. As in Fig- 
ure 1, the dotted line in each panel is the apparent speed for a 
jet with a bulk Lorentz factor of 30 as a function of 6. If all 
sources had the same Ti„, and none had T > 30, the dotted line 
would encompass all of the sources. 

Figure 2b is a plot of Papp versus the maximum observed 
brightness temperature for each source in our real sample. 
The solid line represents sources at the critical angle which 
have Tint = 2 x 10" K, and the same value for 7;„, is used 
in calculating the dotted "envelope." This value of Ti„, was 
chosen so that approximately 75% of the sources would fall 
below and to the right of the solid line; thus matching our 
expectation from the simulation. Figure 2c is a similar plot 
except that it shows sources in their median-low state, what 
we call the "25% Median". This 25% median is the median 
of the lowest half of the brightness temperature observations 
for a given source. In this way, the 25%-median represents 
a typical low brightness state for each source. The ~ 75% 
line for the 25%-median values is !;„, = 3 x 10'" K. A more 
detailed fit of our data to the distribution of simulation points 
in Figure 2a is beyond the scope of this paper, and we note 
that other reasonable choices for the simulation parameters 
described in section |2] will give a roughly similar distribution 
of points and are likely to yield fractions between 60% and 
80% of sources within the critical angle. For a 60% fraction, 
we would have deduced Ti„, = 3 x 10" for the maximum case 
and Tj,,, = 4 x 10'" K for the median-low case. 

It is interesting that plotting f3app versus T„bs in figures 2b 
and 2c produces a trend very similar to that seen in the simu- 
lation, although with more scatter. There is clearly a lack of 
sources with both low observed brightness temperature and 
fast proper motion. The fact that the trend is still clear here 
implies that both plots are described by a narrow range of in- 
trinsic brightness temperature, centered on the value chosen 
for the solid line in each plot. If the scatter in the 25%-median 
and maximum plots is due to intrinsic brightness temperature 
alone, we estimate the real range of intrinsic brightness tem- 
perature in each case extends down to about 50% of the 7},,, 
for the plotted line and up to about twice this value. Note 
that this means the intrinsic brightness temperature ranges for 
the 25%-median and maximum cases do not overlap. Real 
sources appear to be in a distinctly different state when they 
display their maximum brightness temperatures. 

We note that the apparent increase in brightness temper- 
ature between the 25%-median and maximum states cannot 
be attributed to changes in the jet Doppler factors, which we 
have assumed to be constant in time, because the necessary 
large changes in Doppler factor would produce much larger 
changes in the flux densities of the jet cores (oc 6^) which are 
not observed. 
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Fig. 2. — Plots of apparent motion versus observed brightness temperature. 
As described in the text, panel (a) contains 106 sources from the simulation 
illustrated in Figure 1, choosen to simulate the selection characteristics of 
our real data. Panels (b) and (c) plot our real data. Lower limits are indicated 
by arrows and solid dots represent measurements. The solid line represents 
sources observed at the critical angle that have the intrinsic brightness tem- 
perature indicated in the upper left hand comer of the panel. The dashed 
line represents the possible apparent speeds of a F = 30 source with intrinsic 
brightness temperature given by 7],,,. 



The maximum plot contains more lower limits than mea- 
surements, so the characteristic value of 7^,,, should really be 
taken as a lower limit on the characteristic value. Therefore, 
Tj,,, > 2 x 10" K when sources display their maximum bright- 
ness temperature. By equation 1, this requires sources in that 
state to be well out of equipartition, with 77 ^ 10^ or greater 
On the other hand, the characteristic value for the median- 
low state is much closer to the equipartition brightness tem- 
perature defined by Readhead ( 1994), and there are far fewer 
limits in that plot'"', so it seems reasonable to take the char- 

''^ We note that we tried to examine the pure median case, as representative 
of the average state of a source; however, even the pure median case had a 
large number of limits (about half), and we were therefore not able to do a 
robust comparison with the maximum case. The 75% line for the pure median 
case was at 6 x lO'^/f. 
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acteristic median-low intrinsic brightness temperature to be 
Tin, =i 3 X 10'° K. 

Our data indicate that AGN jet cores are indeed near 
equipartition in their median-low state; however they go well 
out of equipartition in their maximum brightness temperature 
state. In their maximum state they have > 10^ times the en- 
ergy in their radiating particles as they have in the magnetic 
fields. Such a large imbalance may be the result of injection 
or acceleration of particles at the base of the jet, resulting in 
a transient state that is in excess of the equilibrium equipar- 
tition brightness temperature. For 77 ^ 10^, equation 5 from 
iReadhead ( 1994) predicts a particle number density which 
has increased by a factor of 400 and a magnetic field strength 
which has decreased by a factor of 15 compared to the equi- 
librium value. We emphase that these values are only approx- 
imate as the precise balance between particle and field energy 
is a sensitive function of brightness temperature (see equation 
1). What is robust, is that the jets in our sample have a narrow 
range of intrinsic brightness temperatures in their median-low 
state, and the brightness temperatures in their maximum state 
are nearly an order of magnitude (or more) larger, indicating 
a much higher particle to field energy ratio. 

It is interesting to note that the narrow range of intrinsic 
brightness temperatures when sources are in their median-low 
state means we can derive Doppler factors for most sources 



which are good to a factor of two or better by assuming 
Tin, = 3 X 10'" K for the 25%-median observed value. Then 5 = 
Tobs(25%-median)l'i X lO'" K. This is Very similar to the equipar- 
tition Dopp ler factor suggested by Sinaal & Gopal-Krishn^ 
("1985) and Readhead ( 1994); however, by using the 25%- 
median observed brightness temperature we get a much more 
reliable estimate for the Doppler factor. Doppler factors de- 
rived from single epoch brightness temperature measurments 
(e.g. Zensus et al. 2002) are likely to have a large degree of 
scatter as it is clear from this work that the intrinsic bright- 
ness temperatures of sources can vary by an order of magni- 
tude or more. Ad ditionally, Doppler factors de rived from out- 
burst events (e.g. lLahteenmaki & Valtaoiall99 9) are likely to 
be systematically overestimated if it is assumed that sources 
are near equipartition during these events. 
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